Heat Bath Approach to Landau Damping and Pomeranchuk Quantum Critical Points 



Johan Nilsson and A. H. Castro Neto 
Department of Physics, Boston University, Boston, Massachusetts 02215 
(Dated: February 2, 2008) 

We study the problem of the damping of collective modes close to a Pomeranchuk quantum 
critical point in a Fermi liquid. In analogy with problems in dissipative open quantum systems, we 
derive the Landau damping of a Fermi liquid by integrating out a macroscopic number of degrees of 
freedom from a generating functional. Being a reformulation of the linearized Boltzmann equation 
this approach reproduces well-known results from the theory of Fermi liquids. We also study the 
Bethe-Salpeter equations within the Landau theory and discuss the implications of these results on 
quantum phase transitions of the Pomeranchuk type and its dynamical exponent, z. We apply our 
results to the electronic nematic instability and find z = 3 in the collisionless limit. 
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I. INTRODUCTION 

Recently there has been interest in the so-called quan- 
tum critical Pomeranchuk instability where the interac- 
tions in a fermion system are strong enough to make the 
Fermi surface soft and ultimately unstable toward a spon- 
taneous deformation^. These instabilities have been pro- 
posed to be behind exotic electronic phases such as elec- 
tronic liquid crystals in cuprate superconductors^* 3 - and 
fractional quantum Hall effecli^, metamagnetic transi- 
tions in electronic systems such as Sr 3 Ru20^, electronic 
systems close to van Hove singularities^, "hidden or- 
der" in URu2Si2 8,9 , a mechanism for generation of spin- 
orbit coupling™, and even phase transitions in interact- 
ing quantum dotsAi. The Pomeranchuk instability can 
take place both on a lattice and in the continuu m 2 ' 12113 
and is driven by attractive quasiparticle-quasiparticle in- 
teractions. These interactions in Fermi liquid theory are 
parameterized by the Landau parameters, F^ a , where s 
(a) stands for the symmetric (anti-symmetric) spin chan- 
nel and m = 0, 1, ... for the angular momentum (we fo- 
cus on the case of the two-dimensional isotropic liquid 
although the three-dimensional case is completely analo- 
gous). Pomeranchuk showed that if some F^ a becomes 
too attractive {F^f < — 1 in two dimensions) at T = 
then the Fermi surface becomes unstable and cannot sus- 
tain collective oscillations. The simplest example of such 
a transition is a ferromagnetic (or Stoner) instability that 
occurs when Fq < — 1, which is characterized by a diverg- 
ing magnetic susceptibility, x (x 1/(1 + Fq)) A similar 
instability with a diverging compressibility, /C, can be 
achieved by pulling on a liquid to put it into a negative 
pressure region when Fq < —1 (JC oc 1/(1+Fg))14. There 
is also the possibility that the instability takes place in a 
higher angular momentum channel. For example, a de- 
formation into an elliptical Fermi surface can be achieved 
by having a strong enough attraction in the d-channcl, 
that is, F£ < —1. The broken symmetry phase in this 
case has been proposed to be an electronic nematic be- 
cause of the nomenclature in the classical theory of liquid 
crystals where an analogous phase exists. The properties 
of the broken symmetry phase has been found to be very 



peculiar with marked non-Fermi liquid properties^. Al- 
though Fermi liquid theory cannot predict what is the 
nature of the broken symmetry phase, it can be used 
to study the approach to the instability from the Fermi 
liquid (quantum disordered) side of the transition (see 
FigQJ. 
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FIG. 1: Phase diagram of a quantum critical Pomeranchuk 
instability of order m. Vertical axis: temperature T; Hori- 
zontal axis: Landau parameter F m . The dashed arrow line 
shows the direction of approach to the quantum critical point 
used in this work. 

One of the main characteristics of a Pomeranchuk 
quantum critical point is that in its proximity the excita- 
tions are not of the simple quasiparticle type but rather 
become heavily damped because of the soft nature of the 
Fermi surface. This is so because the Fermi surface allows 
for a diverging number of low energy particle-hole pairs 
to accompany an excited quasiparticle. Presumably, a 
similar situation holds for most of the single-particle ex- 
citations in the broken symmetry phased although the 
Fermi liquid theory is not able to access that phase di- 
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rectly. 

A natural way to study the collective modes of an in- 
teracting electron fluid is via Fermi surface bosonization. 
The method of bosonization has been very successful in 
providing a framework for describing the breakdown of 
the normal Fermi liquid in one-dimensional (ID) systems, 
where instead the Luttinger liquid emerges as the "stan- 
dard" modeli^. In Luttinger liquids there are no well 
defined quasiparticles and the electron breaks down into 
charge and spin collective modes. It is now well under- 
stood that the anomalous behavior in one dimension is 
mainly due to the restricted phase space. Indeed, there 
exists Pomeranchuk instabilities also in ID which lead 
to a breakdown of the Luttinger liquid. What happens 
is that either the charge- or spin- velocity vanishes at 
the transition. Since the charge(spin) susceptibility goes 
as the inverse of the charge(spin) velocitjsi^ the signa- 
ture of the instability is a diverging susceptibility, just as 
in the higher dimensional case. Thus, the Pomeranchuk 
instability is a general concept, independent of the di- 
mensionality of the system under interest. The higher 
dimensional bosonization approach was developed many 
years agoiLiLigiigiSSiSi and reproduces the known results 
from the theory of Fermi liquids, especially the physics 
of collective modes and particle-hole excitations which 
are described in terms of boson coherent states. The ba- 
sic reason for this is that both theories describes fermion 
systems dominated by forward scattering correctly, and 
this is what the Landau Fermi liquid is all about^. 

For problems such as those involving classical limits it 
is a good idea to use a description in terms of bosons since 
they can have classical analogs whereas fermions cannot. 
In particular, when one is studying phase transitions in 
the language of a Landau-Ginzburg expansion the order 
parameter is a classical construction which in the quan- 
tum language should be described by bosons. One can 
generally trade a fermion system with density-density in- 
teractions to one involving bosons by using a Hubbard- 
Stratanovich transformation and then integrate out the 
fermions from the path integral that describes the sys- 
tem. This is the standard device to study quantum criti- 
cal theories in Fermi systems22*2£i2£. In such approaches 
to quantum critical phenomena dissipative terms are gen- 
erated in the order parameter effective action, the so- 
called Landau damping. However, there have been some 
arguments about whether it is viable to integrate out all 
of the fermions from the system in this way, or if some 
properties of the system is lost due to the gapless nature 
of the electrons in a Fermi liquid. In particular there are 
particle-hole- like excitations with arbitrarily low energies 
which must be treated with care since they may control 
the critical behavior of a Fermi system. Nevertheless, 
the integration of high energy degrees of freedom as one 
usually does in renormalization group (RG) approaches 
should give sensible results. Obviously, the particle-hole 
continuum type excitations do not disappear in a faith- 
ful bosonic reformulation of the system, although their 
nature might be disguised. When some of these modes 



are in resonance with the collective excitations the result 
is Landau damping. This occurs only for modes with a 
velocity smaller than the Fermi velocity. 

In this paper we show that the Landau damping prob- 
lem can be recast in terms of a trace over bosonic modes. 
In particular the modes associated with the high angular 
momentum oscillations of the Fermi surface are traced 
out since the collective excitations we are interested are 
concentrated in the low angular momentum modes. Fur- 
thermore, using the Bethe-Salpeter equation in the con- 
text of the collision integral, we can study the lifetime of 
quasiparticles and collective modes in a Fermi liquid close 
to a Pomeranchuk quantum critical point. Moreover, we 
discuss the dynamical critical exponent of the same quan- 
tum critical point. The paper is organized as follows. In 
Section ^ we show how Landau damping as described 
above can be reproduced within a simple purely bosonic 
reformulation of the Fermi liquid theory. Here the damp- 
ing comes about from the coupling to a bath of harmonic 
oscillators similarly to what happens in dissipative quan- 
tum mechanics. In Section IIIII we study the collision 
integral in the language of Fermi liquid theory, focusing 
on the behavior near the critical point. This turns out 
to be important since collisions sets the lifetime of the 
quasiparticles, and as one approaches the critical point 
the lifetime decays rapidly. We conclude in Section IIVI 
with a discussion of the implications of these phenomena 
on the theory close to the critical point. In an Appendix 
we present some mathematical details. 



II. LANDAU DAMPING FROM A BOSON 
BATH 

The phenomenon of Landau damping is conventionally 
thought of as the damping of a collective mode due to 'a 
coherent interaction of the collective mode with those 
particles which "surf-ride" on the crests of the running 
wave. — In this section we show explicitly how Landau 
damping can also be viewed as the dissipation of a mode 
due to its coupling to a set of independent oscillators, 
i.e. a "boson bath" . Here the origin of these bosons 
are particle-hole excitations of the fermion system itself. 
This approach is similar in spirit to theories of dissipation 
in quantum mechanics 

Let us take as our starting point the Landau kinetic 
equation that describes the evolution of the quasiparti- 
cle distribution function in the spin symmetric channel, 
np{r, fp.: 

— ^ h V p ep(f, t) ■ V r np(r, t) 

- V r ep-(r, t) ■ Wpnp(r, t) = I[np>] , (1) 

where ep(r,t) is the quasiparticle energy and /[/ip»] is 
the collision integral (p is the quasiparticle momentum, 
r is the position in real space, and t is time)2S. Since the 
physics goes on close to the Fermi surface it is convenient 
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to change variables to vp{f, t) which we define by writing 
the deviation from the ground state values of quasiparti- 

cle distribution as 5np(f,t) = — yjj^)V p {r,t). Physically 
v p (r,t) measure the local deformation of the Fermi sur- 
face. Upon a Fourier transform in space the linearized 
kinetic equation for v p (r , t) becomes 



dvpjg, t) 
dt 



+ ivp-q [up(q, t) + 5ep(q, t)] = l'[vp>], (2) 



where vp is the Fermi velocity at the momentum p. Again 
the notation 



dtp ^ fpp fififp j 



(3) 



is standard and describes the interaction (parametrized 
by fpp') among excited quasiparticles. For an isotropic 
liquid in two dimensions, vp is parallel to p, and since the 
magnitude of p is nearly equal to the Fermi momentum 
Pf we can label p by an angle 9 P (measured from the 
direction of q). Then, the interaction fpp only depends 
on the angle between the momenta and is convenient to 
expand in Fourier components because of the circular 
symmetry: 



f(0,0')^J2 m fme^ e - e '\ 



(4) 
(5) 



where f m — f- m - The sum over m should be cut-off at 
some integer of order Ep/A = N\ ^> 1, where E F is the 
Fermi energy and A is some energy scale below which 
the Landau theory is valid. Projecting onto the Fourier 
components and setting the collision integral to zero for 
the time being we arrive at 



du m (q,t) iv F q 



Of 



(1 + F m -i)v m -i(q,t) 
+ {l + F m+1 )v m+l {q,t) =0. (6) 



Here vf is the Fermi velocity and the F n = N(0)f n are 
the interaction parameters normalized with the density of 
states at the Fermi energy, iV(0). Note that these equa- 
tions are linear in the time derivatives and couple time 
derivatives of odd components to the even components 
and vice versa. Thus the equations respect time reversal 
invariance implying non-dissipative dynamics. 

To make connection with the more familiar simple os- 
cillator we eliminate the odd components in favor of a 
second time derivative 30 . With the notations: 



a m (q,t) 

B m 



yl + F m v m {q,t) , 



(1 + F m )(2 + F m _i + F m+1 ). 



(7) 
(8) 



y/1 + F m _i(l + F m )y/1 + F rn+1 , (9) 



the equations of motion become: 
d 2 a m (q,t) 



dt 2 



v F qY A 



+ B m -ia m -2(q, t) + B {q,t) =0. (10) 

These equations of motion are also generated from Euler- 
Lagrange equations and the classical Lagrangian density: 



^ a 2m(g^) da 2m (q,t) 
dt dt 

- ( v Fq) 2 ^2a2 m (q,t)V mn a 2n (q,t), (11) 



where 

2m-l + <5m+l,n-E>2 m +l . (12) 

Notice that (|llf> describes a set of coupled harmonic os- 
cillators coupled by " springs" with spring constants given 
by Vm.n- This result shows the clear collective nature of 
the Fermi surface modes. 

We can also encode the evolution of the system with 
the methods of path integrals 31 in the real time partition 
function: 



Z 



Da* Da e 



is 



(13) 



where we also perform a Fourier transform in time 



S = 



q,m,n 



du>a* 2m {q,Lo)V mn a 2 n{q,^), (14) 



with, 



Vrr, 



(v F q) 2 [s 2 ^ 
uj/(v F q) . 



(15) 
(16) 



It is clear from these results that the propagator of the 
collective modes is essentially given by the inverse of the 
coupling matrix V m , n , that is, V^ n . Therefore, the eigen- 

modes of V are the collective excitations of the system 
and its eigenenergies determine the frequency of oscilla- 
tions of these modes. 

We can gain further insight into the problem by split- 
ting the action 114|) into two decoupled sectors (S = 
S + + S~) by using the symmetric and antisymmetric 
combinations 



1 

71 



(17) 



defined for m > 0. These corresponds to the longitudi- 
nal (+) and transverse (— ) modes in the sense of having 
even/odd parity with respect to the line defined by q. 
Notice that the transverse modes have no density (i.e. 
m = 0) component. 
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Formally, there can be an infinite number of Landau 
parameters F m that parameterize the interactions in a 
Fermi liquid. In practice, however, only a few Landau 
parameters are taken into account* 2 ^ In the case of the 
Pomeranchuk instability we can concentrate on a single 
Landau parameter, that is, the most singular, since all 
the other parameters only provide trivial renormaliza- 
tions of the Fermi liquid properties. Without lack of 
generality let us assume that there are M Landau pa- 
rameters in the description of a quantum liquid, that is, 
the F m are zero for m > M + 1. In this case the matrix 
V m ,n hi JTSJl for m,n > M + 1 (we call it the high-M 
block) can be written in the longitudinal and transverse 
sectors in terms of blocks of the form: 



(s 2 - 







V 



A 



7 



(18) 



The matrix (|18fl (assumed to have size N x N) have or- 
thonormal eigenvectors of the form: 



<p(k) = \j N 2 +l (sin(fc), sin(2fc), sin(3fc), . . . , sin(7Vfc)) T ; 

(19) 

with respective eigenvalues A(s, k) — s 2 — (1 + cos(fc))/2 
(the allowed values for the fc's are: k = ixnj (JV+1) where 
n = 1...N). The block (JTHJ is coupled to the mode M 
via two off-diagonal terms in V m ,n- Schematically we 
have: 



\ 



low-M block 

• ■ ■ B 2 m+i 
■■■ 







••• 



B2M+1 



high-M block 



V 



(20) 



Using the diagonalized form of the block matrix it is not 
hard to integrate out the high-M modes. The end result 
is to change 



A 



M 



A 



M 



8 A 



M ■ 



in the low-M block of V mn , where 



and 



G(s) 



SA M = 2B 2M+1 G(s) 

sin 2 (A;) 



N ■ 



- (l + cos(fc))/2' 



(21) 



(22) 



(23) 



the limit of N — > 00 the poles become dense resulting in 
a branch cut. By taking N — > 00, G(s) can be expressed 
as an integral: 



G(s) = 



71 dk 



sin 2 (fc) 



% 7T 2,s 2 — 1 — cos(fc) 



(24) 



Taking s to lie away from the branch cut one can turn 
the integral into a contour integral around the unit circle 
centered at the origin. The resulting integral is easily 
evaluated for \s\ > 1 by the residue theorem and can be 
analytically continued to |s| < 1 with the result: 



G(s) = -2 + 4s 2 - 4 v / s 4 -s 2 . 



So that for real s and \s\ > 1 



G{s) =-2 + 4s 2 - A\sWs 2 - 1, 
while for real s and \s\ < 1 



G(s ± 



-2 + As 2 =F 4is\/l - s 2 . 



(25) 



(26) 



(27) 



Hence, the final result of integrating out the high angu- 
lar momentum modes is the introduction of a shift in the 
Landau parameters with the introduction of an imagi- 
nary part when \s\ = \u>\/(qvF) < 1, which is nothing 
but the Landau damping. By only keeping the low-M 
modes the time reversal invariance has been broken. 

Since what we have done is a reformulation of the Lan- 
dau theory it should give the same result as the conven- 
tional approach. A simple check amounts to taking only 
Fq ^ and integrate out all modes except the zeroth 
one in the longitudinal sector, that is, M = 0. The well- 
known zero sound mode should then come out as the zero 
of the only remaining matrix element (which is equiva- 
lent to having a pole in the propagator). Explicitly we 
obtain: 



A + 2d A 



= 0. 



(28) 



Substituting the expression for A , Bi and G(s) we get 
an undamped solution for _Fo > with 



1 



V(l + ^o) 2 



r 



(29) 



which is the dispersion of the zero mode in two dimen- 
sions, see e.g. Ref . [l9| . 



III. LIFETIME EFFECTS CLOSE TO A 
POMERANCHUCK QUANTUM CRITICAL 
POINT 



From this result we see that G(s) as a function of complex 
s is analytic in the upper and lower half plane. It has a 
set of poles on the real axis between s = — 1 and s = 1. In 



In this section we investigate the effect of the prox- 
imity of a Pomeranchuk quantum critical point on the 
single particle and collective modes of a Fermi liquid. 
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Using the collision integral formalism we study the quasi- 
particle lifetime via the Bethe-Salpeter equation close to 
the quantum critical point. We show that the critical 
point associated with the quantum fluctuations change 
substantially the quasiparticlc lifetime without affecting 
the damping of the collective modes at low angular mo- 
mentum. 



A. Single particle states 

We follow the usual procedure and study the lowest 
order process, i.e. quasiparticle-quasiparticle scattering. 
Then the expression for the collision integral becomes^: 



clarity we then get 

; = ra/""££ |<34|t|12> ' 2 

q p> 

8(vFqcos(6 p ) — uj)S(vf qcos(9 p >) — u) 

[1 - n(e - u))]n{e') [1 - n(e' + u)] . (33) 

This form is convenient when e > 0, for e < one can 
instead exchange n <-> 1 — n everywhere. The end result 
is even in e. 

To proceed we need an expression for the matrix ele- 
ment. This one can get from the solution of the Bethe- 
Salpeter equation, which describes repeated scattering of 
a quasiparticle-quasihole pair: 



/K-] = 27r^]T|<34|t|12>| 2 

P2&2 P3 CT 3 

%i+P2,P3+?4^1+CT2,tT 3 +(T4^( e l + £ 2 — ^3 — e i) 

n 3 n 4 (l - n 2 )(l - ni) - nin 2 (l - n 3 )(l - n 4 ) , (30) 



which describe scattering processes in and out of the 
state p which is labeled 1. The spin information can 
be encoded into the the definition of the matrix element 
| < t > | 2 as described in Ref. j2^, so we drop the spin 
information from now on (or assume a spin-independent 
scattering amplitude). 

To make progress we linearize the collision integral in 
the deviations 6rii from (local) equilibrium n°. More- 
over, we employ the relaxation time approximation which 
amounts to keeping only the term proportional to 8n\: 



I 



1 



-Snp, 



(31) 



where 



2^EE l<34|t|12>| 2 



VI P3.P4 



P1+P2.P3+P4 



sji^i + e 2 - e 3 - e 4 ) 



n§n2(l-n§)+n§(l-ng)(l-n2) 



(32) 



One can do better than the relaxation time approxima- 
tion but it is not necessary for our purposes^ . 

Now we introduce the momentum transfer q, and the 
energy transfer to via the identity 5(e± + e 2 — e 3 — e 4 ) = 
f du6(ei — e 3 — uj)5(e 2 — e 4 + uS). We also specialize to 
nearly forward scattering. This is an approximation in 
three dimensions where it amounts to neglecting out of 
plane scattering. In two dimensions it is essentially ex- 
act because of momentum conservation and the reduced 
phase space. This means that we restrict the sum over 
q to q < q c , where q c <C pf is some cut-off momentum. 
Dropping the superscripts on the occupation numbers for 



Jpp"~7V~ /^^ t p"p>(q, LJ + ir]). 



p"^p' 



oj + irj — q ■ vp/i 



The formalism for solving the Bethe-Salpeter equation in 
three dimensions can be found in e.g. Ref. |3^ . Here we 
outline the procedure in two dimensions. 
Firstly one expands t in Fourier components 

t&%u> + iri)= J2 t ^m'(q,^ + iv)e iimep - m '^' ) , (34) 



where the angles 9 are measured with respect to q. A 
matrix equation for t m , m ' can be derived from the or- 
thogonality between the Fourier components 

(1 ~\~ F m ^t m m ' — m' fm ~1~ Fjji ^ K(n 17l)t nm > 7 (35) 

a 

where 



K{n) 



d6 



o 2ns — cos(0) 



and 



uj + %r\ 
VFq 



(36) 



(37) 



Explicit expressions for Kin) can be found in Ref. |19j . 
Given the values of {/„} one can solve <|55|) to get t m>m >. 

Close to a Pomeranchuk quantum critical point when 
one of the Landau parameters reaches a critical value we 
can keep only the relevant Landau parameter (see Fig^). 
In this paper we study the "nematic" critical point and 
choose to keep only F 2 — F_ 2 nonzero in which case 
it is straightforward to solve for t. The only nonzero 
components are 



_ 1 + F 2 - F 2 K(0) 

E2.2 — 1-2. -2 — Y) -> 2 • 

_ _ F 2 K{A) 

52,-2 — 1-2,2 — Y) J2 , 



(38) 
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where 



D = [1 + F 2 - F 2 K(0)} 2 - [F 2 K (4)] 5 



(39) 



Now that we have an explicit expression for < t > we 
can go back to the expression for the lifetime in l|33|l , turn 
the sums into integrals and perform the angular integrals 
with the result 



N(0)L 2 



(2irv F ) 2 [1 - n(e)] 



J duj [1 — n(e 



x I de'n(e')[l-n(e' + u)} F ^\t(s,q)\ 2 , 

J\u\/v P 1 



(40) 



where 

l^,9)| 2 



16 



■[l-85 2 (l-S 2 )] \t 2 , 2 +t 2 ,- 2 \ 



+512 S 4 (l- S 2 )(|f, 2 , 2 | 2 + |i 2 ,_ 2 | 2 ). (41) 



The first term in l|41|l gives the dominating contribution 
(from small s) so we only keep it in what follows. Notice 
that due to the integration limits on q in l|4UI) we have 
s < 1 and the singularity near s = 1 is removed by the 
behavior of \t 2 _ 2 + t 2 ~ 2 \ 2 . 

Usually, i.e. for not too strong attractive interaction, 
it is a good approximation to take t to only depend on s. 
Close to the critical point however, one must also include 
a g-dependence in t, explicitly we take 



1 + F 2 = S + n(q/q c 



(42) 



where k oc d 2 F 2 /dq 2 . Now we can trade the integral over 
q to one over s. The integral is then dominated by the 
small- s contribution and it is possible to approximate the 
integral by a scaling argument. For small s we have: 



^2,2 + t 2 . 



-2 — 



h 



F 2 -F 2 [K(0)+K(4)} ' 

5 + n(q/q c ) 2 + 2iF 2 s ' 

5s 2 + {uj/uj c ) 2 -2is 3 ' 
Using l|43[) the g-integral in (|4(J|) can be written as: 

,,3 



16|/ 2 | 2 



ds 



M/<w«. l$s 2 + (^K) 2 ] 2 + 4s 6 ' 



(43) 



(44) 



For 5 — * (that is, at the quantum critical point) we 
scale s 3 oc uj 2 so that the integral becomes oc cl>~ 4 / 3 . 
Performing the remaining integrals at T = one get 



1 



— oc e 

T 



2/3 



(45) 



leading to a very anomalous energy dependence of the 
quasiparticle lifetim e) 2 ! 12 . For finite S (away from the 
quantum critical point, inside of the quantum disordered 



regime) we should recover the standard Fermi liquid re- 
sult. By scaling socu, we get: 

- oc e 2 \n(u c /e) 2 + const x e 2 , (46) 
r 

which is the Fermi liquid resulli 1 ^. The above results can 
be summarized by the scaling form of the quasiparticle 
lifetime that can be obtained using the contour integral 
trick of Pethick 3 ^ that we give in the appendix. The 
general result for the quasiparticle lifetime is: 



. (47) 



where F(£) is a scaling function such that F(£) w 1 for 
£ 3> 1 and w — £ 2 / 3 log(£) for ( < 1, reproducing 

the results of eqs. Ij45(l and (|46|l . A similar crossover 
formula (for the imaginary part of the self-energy) in 
the case of a ferromagnetic quantum critical point was 
found by Chubukov et al£&2^. Notice that our results 
can be generalized for any kind of Pomeranchuk quan- 
tum critical point given the critical Landau parameter. 
The result of this lowest order approximation is thus that 
close to the "nematic" critical point in 2D the lifetime of 
the quasiparticle dominates over the energy and the con- 
cept of the usual quasiparticle is not well-defined. To see 
the real fate of the quasiparticle one must include higher 
order scattering processes or use a more microscopic the- 
ory. Nevertheless, the calculation in this section signals 
the breakdown of normal Fermi liquid behavior close to 
a Pomeranchuk quantum critical point. 

Studies by other groups of similar models also find a 
breakdown of Fermi liquid theory at the critical point. 
The e 2 / 3 behavior in two dimensions has also been found 
before in related systems, see e.g. Ref. 0E 



B. Collective modes 

To calculate the relaxation time for collective modes 
one must keep in mind that inside the linearized collision 
integral there are terms of the form 5np — 8np+^ where 
the second term is neglected in the usual relaxation time 
approximation. The existence of the second term makes 
collisions less effective in the damping of collective modes. 
More explicitly, linearizing the collision integral (|30[l for 
the v's in <|2J) we get: 



I'M = - 



2n 



£]T|<34|£|12>| 2 



P2 P3.P4 



P1+P2-P3+P4 



?4 5(ei + e 2 - €3 - e 4 ) 



n 2 (l - n 3 )(l - n 4 ) v\ - 1/3 + v 2 - v± 



(48) 



Expanding the v J s in Fourier components as in Section 
Inland projecting the collision integral we get 



1 



1 



(49) 
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where 1/r' is given by (|40|l with the substitution J dq — > 
/ q 2 /2p 2 F . To get this result we have neglected the 
term — since this term is typically of the same size 
as the first (i.e. ^1 — ^3) and just gives a different prefactor 
in the calculation of 1/r of order 1. More important is 
the expansion of cos[m(# + q/pF)] in powers of mq/pF, 
which is only valid for m <C Pf/Qc ~ N\. The terms 
with odd powers of q vanish upon performing the angle 
averages. 

For these modes the simple scaling evaluation of the 
integrals gives: 



for the transverse 



mode and 



oc 



/ m 
\N~ A J 



for (5 — > and 



(50) 



(51) 



for finite S. From this we can conclude that the collision- 
less approximation is good for self-driven modes (ui ~ e) 
when (m/Nh) 2 <C 1. Note that this implies that it makes 
sense to talk about (e.g.) the M = 2 mode also in the 
critical regime. The high-Af modes are more sensitive 
to low-<7 scattering. For these it is not possible to re- 
late cos[m(# + q/pp)] in any simple way by cos(m0) and 
hence the single particle result above should be a good 
approximation. This means that the high-M modes are 
heavily damped in the critical regime. For intermediate 
M we expect a crossover between the two limiting cases. 
Also, note that in all cases (single-particle or collective) 
there is a prefactor of l/(kpVF ) ~ 1/cf an d the correct 
powers of lu c to make the dimensions correct. 



IV. QUANTUM CRITICAL REGIME - 
DISCUSSION 

Recently Yangi^ has argued in favor of a dynamical 
critical exponent of z = 2 for the Pomeranchuk quantum 
critical point in 2D. This is in disagreement with e.g. the 
results of Ref. Q . The claim that the other non-critical 
modes just renormalize the couplings in the critical the- 
ory is not obvious. In particular, the low-energy particle- 
hole fermion excitations which are reproduced within the 
bosonic formulation as linear combinations of the high 
angular components are neglected in Ref. 01 ■ The cal- 
culation in Section ITT1 illustrates the effect of these modes 
in the collisionless regime, and in a simpler setting than 
that of multidimensional bosonization which requires a 
complicated Bogoliubov transformation to achieve simi- 
lar results^ 

Taking only Fq, F\ and F2 ^ we can get the equa- 
tions for the transverse and longitudinal modes from the 
result of Section [HI If we integrate out all the modes 
that we are not interested in (i.e. M > 3), the resulting 
equations for the modes are: 



(s% - Ao)/2 -B 1 



Ao - SA-z 



(53) 



for the longitudinal (+) mode. From these equations 
one can extract the leading scaling behavior to be s 2 _ ~ 
[5 + K(q/q c ) 2 ]/4 and s+ ~ —i[5 + n{q/q c ) 2 /2], assuming 
that only the <i-channel goes critical. So the transverse 
mode has a dynamical critical exponent of z = 2 a nd the 
longitudinal z = 3 in agreement with Ref. |2l37j |. The 
longitudinal mode, being slower, dominates the critical 
behavior. 

As we found in Section ITTT1 the collisionless behavior 
is only a good approximation for low-Af modes in the 
critical regime. The question that arises is what would 
be the effect of collisions on the modes that participate 
in the Landau damping. The coupling of the M = 2 
mode to the rest of the system still results in damped 
motion. Generally one can integrate out the modes with 
M > 3 with the same result as above in (|52|l and l|53[) . 
The difference is that we no longer have an explicit ex- 
pression for SA2 since the approximations of the previ- 
ous section break down. However, we still must have 
SA2 = 2B 2 U(q,uj), where U(q,u>) is a function such that 
[/(<?, 0) = U(s = 0) = -2. The calculation of U(q,u) is 
beyond of the scope of this work. Nevertheless, if U(q, u>) 
is known, the longitudinal mode equation becomes: 



s 2 -i 

S+ 8 





-i-Y 


X 









U(q,Lu)-U(q,Q) =0. (54) 



A 2 - 8A 2 



0. 



(52) 



If we assume that U(q,ui) — U(q,0) ~ —iuj a /q in the 
lj — > limit, the dynamical critical exponent of the mode 
is z — 3/(2 — a). This result allows us to put bonds 
on the value of the dynamical exponent since we expect 
2/3 < a < 1 in such a way that 9/4 < z < 3. One 
can speculate that the modes that are responsible for the 
damping of the critical mode, having very low energy, are 
weakly damped since there is not much for them to decay 
into because of the small phase space. This would hint to- 
wards z — 3 even when collisions are taken into account. 
Another argument in favor of z — 3 is that the collision- 
less limit is a good approximation when the mean free 
path is much longer than the wavelength, which trans- 
lates into 1/r <C qvF^ Taking 1/r ~ tv 2 / 3 and to ~ q 3 
the inequality is clearly satisfied as q — > indicating the 
validity of the collisionless limit. 

A questionable thing with the whole approach close 
to the critical point is whether the calculations are self- 
consistent. The simplest particle-particle scattering pro- 
cess gave a lifetime that is proportional to e 2 / 3 . The 
same result can be obtained from a 1-loop calculation of 
the imaginary part of the fermion self energy, which de- 
scribes similar physics. Taking this seriously would inval- 
idate the quasiparticle picture, and hence make the whole 
ground for the calculation questionable. Notice, how- 
ever, that although the quasiparticle description seems 
to break down, the low angular momentum collective 
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spectrum remains essentially unchanged (except by the 
damping of the collective modes given by the subleading 
contributions (|50[1 and (|51|) L 

One further circumstantial evidence for a breakdown of 
the quasiparticle near the quantum critical point comes 
from the theory of multi-dimensional bosonization 20 . 
The calculations are tractable for an isotropic interaction, 
the connection to the notations in this paper is Ng = Fq. 
The Bogoliubov transformation used to diagonalize the 
system connects different patches i and I of the Fermi 
surface according to the matrices: 



Mi 



Su + 



1 Ng 



Mi = - 



N 1 + Ng Sl - Si 
1 Ng y/SiSi 
Nl + Ngsi+Si 



(55) 



where s, = cos(0j). Away from the critical point the sec- 
ond term in M. represent the dressing of the particle by 
interactions, and M is always small. Close to the critical 
point however, the factor Ng/(1 + Ng) = F /(l + F ) 
blows up when Fq — ► — 1 and a broad dressing cloud 
dominates the excitations, also N is no longer negligi- 
ble and is an important part of the excitation. More- 
over the expression for the quasiparticle residue Z = 

exp{- 2 } shows that it vanishes exponentially 

as the transition is approached. The calculation for the 
d-channel Pomeranchuk transition is much more compli- 
cated because the kernel of the integral equation is no 
longer a product kernel. We expect similar conclusions 
however. Because the Fermi liquid theory breaks down 
close to the critical point a calculation within the Fermi 
liquid theory itself doesn't make much sense there, except 
for indicating its own breakdown. The correct descrip- 
tion is not in terms of the standard Fermi liquid theory 
and it would be very interesting to know what it is. 

In summary, we have presented an alternative ap- 
proach for the study of Landau damping in Fermi liquids 
based on the trace over high angular momentum modes. 
We find that this approach reproduces well-known re- 
sults in the literature and allows for the study of the 
quasiparticles and collective modes close to a Pomer- 
anchuk quantum critical point. Although our approach 
is general and can be used to study any instability of 
the Pomeranchuk type, we focused on the problem of the 
electronic nematic that has been discussed in the litera- 
ture recentr y 2 i 13 i 38 i 39 i 40 . We found that the quasiparticle 
lifetime becomes anomalously short at the quantum crit- 
ical point leading to a breakdown of the Fermi liquid 
description. 
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APPENDIX: CONTOUR INTEGRAL 

The formula for the inverse lifetime can be written (at 
T = and e > for simplicity) 



1 



1-KVpkp 



dujuj 



-c<+iV 2 (0) / 1 -|i( S , 9 )| 2 



1 



dcouj(-cf + I s ). (A.l) 



To get this we have added and subtracted the contribu- 
tion of the N 2 (0) times the s-integral down to zero, this 
is approximately 



c< = 16 



\/v F qc 



ds 



[ds 2 + {w/uj c ) 2 ] 2 + 4.s 6 



1 ~3 

16 / dx - 

Jo [Sxi + ^ + ix^y 



, (A.2) 



which to leading order is a constant (for finite n), note 
that we set F2 = — 1 everywhere where it is not dangerous 
to do so. 

Now we can follow Pethick 3 ^ and turn the integral in 
I s into two contour integrals, the difference being that in 
our case the two contours have different orientation 



2N(0) 



1 



s 2 v / T — ^ 2 l + 4s 2 (.s 2 -l) 



(*2,2 + ^2,-2)) 

(A.3) 



using the fact that the integral is dominated by small 
values of s one can rewrite it as 

^-^r ^-fr./Li'+h- <a - 4) 

after a deformation of the contour to the imaginary axis. 
Scaling z — 5x/6 we get 



= -16(f) 9 ? 



dx- 



1 



x" + 3x 2 -4Y } 



(A.5) 



where Y = (3/5) 3 (ui/uj c ) 2 . Finally performing the re- 
maining integral we obtain 



1 



1 



1 



+ 2e 2 (^-^V J dxlog 



3x + 2(1 -20 



x 3 - 3x + 2 



(A.6) 
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which is of the desired scaling form with £ = 
(3/8) 3 (e/u; c ) 2 . 
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